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^^ Abstract. We experimentally measure the monochromatic transmission matrix (TM) of an 

optical multiple scattering medium using a spatial light modulator together with a phase- 

1"^ shifting interferometry measurement method. The TM contains all information needed to 

Q shape the scattered output field at will or to detect an image through the medium. We 

confront theory and experiment for these applications and we study the effect of noise on 
the reconstruction method. We also extracted from the TM informations about the statistical 
properties of the medium and the light transport whitin it. In particular, we are able to isolate 



Q the contributions of the Memory Effect (ME) and measure its attenuation length. 
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Introduction 

For a long time, wave propagation was only studied in the simple cases of homogeneous 
media or simple scattering were the trajectories of waves can be exactly determined [IJ. 
Nevertheless, the effect of the multiple elastic scattering on propagation of waves is complex 
but deterministic, and information conveyed inside such media is shuffled but conserved ^. 
Within the last decades, the study of complex and random media has been an interdisciplinary 
subject of interest, ranging from solid state physics to optics via acoustics, electromagnetism, 
telecommunication and seismology [[3l|4l[5l[6l|71. 

Recent works in optics have shown that it is possible to achieve wavefront shaping beyond 
a complex media despite multiple scattering [8l HI and even to take advantage of the 
randomness of such media to increase resolution fTOl| by reducing the size of the focal 
spots. These works amount to performing phase conjugation. Nevertheless, in acoustics [11 J, 
seismic [fT2l . electromagnetics or telecommunication lfT3ll . other operations are done to 
transmit information through a random propagation systems. These techniques assume that 
the propagation between input and output is characterized in amplitude and phase, which is 
not trivial in optics. 

We recently introduced a method to experimentally measure the monochromatic transmission 
matrix (TM) of a complex media [il4il and presented an application to image transmission 
through such media ifTSl . In the present paper, we will develop the TM model for optical 
systems, in particular for multiple scattering media and detail the experimental acquisition 
method. We will show that the TM can be exploited for different applications. The most 
straightforward ones are focusing and image detection using well-established operators used 
for inverse problems in various domains (telecommunication lfT3]| . seismology [[121 , optical 
tomography [[I6l[l7l|, acoustics [[l8l[T9]] and electromagnetics [|20l ). We will confront theory 
and experiment for applications and study the statistical properties of the TM. We will in 
particular study the effect of the ballistic contributions on the statistics of the TM and its 
consequences on image transmission. We will finally study how the short range correlation of 
the Memory Effect [21 J (ME) affects the TM and how to quantify its correlation length with 
the knowledge of the TM. 

1. Matrix Model and Acquisition 

1.1. Modeling Transmission Through a Linear Optical System 

In free space, a plane wave propagates without beeing modified and its k- vector is conserved 
since plane waves are the eigenmodes of the free space propagation. In contrast, a plane wave 
illuminating a multiple scattering sample gives rise to a seemingly random output field. This 
output field is different for a different input k- vector. In order to use a scattering medium as 
an optical tool, one has to characterize the relation between input and output free modes. At a 
given wavelength, we will model such a system with a matrix linking the k- vector components 
of the output field to those of the input field. We will justify the validity and the interest of 
such a model. 
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For any linear propagation medium, the propagation of an optical wave is entirely 
characterized by its Green function I122II . If the sources are linearly polarized and the 
observation is made on the same polarization, one can use a scalar model. The scalar Green 
function G(r, r', t, f) describes the influence of the optical field at position r' at time f on 
the optical field at position r at time t. For a propagation in a medium stationary over time, 
the time dependence of the Green function is governed hy t — t\ For a surface Sq^c containing 
all the sources, the optical field on r at t reads : 

P P POO 

E{r,t)= G{r,r',t-t')E{r',t')dt'd^r' (1) 

This expression can be written in the spectral domain : 



E{r,uj)= II G{r,r\uj)E{r',uj)d^r' 



(2) 



Where E{r, uj) (resp. G(r , r\ uj)) is the temporal Fourier transform of E{r' ^ i) (resp. 
G{r,r\t) ) at the pulsation uj. We want to characterize the relation between the optical 
field on the surface S^^c containing the sources and the optical field on the elements of an 
observation surface Sohs- Experimentally, sources and receptors have a finite size. We note 
^m * = Lout E{r)dr the average optical field on the m^^ receptor of surface S"^^ C S^obs- In 
a similar way, we note E^^ = L^^ E{r)dr the average optical field on the n^^ controlled area 
of surface S^^ C S'src- We write A^ and M respectively the number of sources and receptors. 

Therefore, we can define the mesoscopic TM of an optical system for a given wavelength 
as the matrix K of the complex coefficients kjnn connecting the optical field (in amplitude and 
phase) at the m^^ output element to the one at the n^^ input element element. Thus, we have : 

N 

^m = 2^ kmnE'^ (3) 

n 

In essence, the TM gives the relationship between input and output pixels at a given 
frequency, notwithstanding the complexity of the propagation for a stationary medium. Along 
this paper, we will perform the Singular Value Decomposition (SVD) of TMs to study their 
transmission properties. The SVD consists in writting : 

K = [/Syt (4) 

Where ^ denotes the transpose conjugate. The SVD decomposes the system in 
independent transmission channels characterized by their transmission coefficients and by 
their corresponding input and output modes. ]/ is a unitary change of basis matrix linking 
input freemodes with transmission channel input modes of the system. S is a diagonal matrix 
containing real and positive values called singular values of K and noted A^ . These values 
are the square root of the energy transmission values of the transmission channels. U is 
the unitary change of basis matrix between transmission channels output modes and output 
freemodes. For practical purposes, we will study the normalized singular values defined by : 
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V Z^k\ 

The S VD is a powerful tool that gives access to the statistical distribution of the incident 
energy injected thought the multiple transmission channels of the system. Within this matrix 
model, the number of channels of the observed system is min{N, M). Nevertheless, this 
number does not quantify the information that could simultaneously be conveyed in the 
medium. 

Quantifying information is of high interest and this work has been done in acoustics Il23l 
[24l for Time Reversal focusing through a multiple scattering medium. In those theories, 
the number of spatial and temporal Degrees of Freedom (A^dof) also called Number of 
Information Grains is introduced. These degrees of freedom represent as many independent 
information grains can be conveyed through the system. For a monochromatic experiment 
there are only spatial degrees of freedom and A^dof is given by the number of uncorrelated 
element used to recover information [2J. This number will be different on the type of 
experiment considered, focusing or image detection. A focusing experiment consists in 
finding the optimal input wavefront to send to maximize energy at a desired position (or 
a combination of positions) on the output of the sample. An image detection experiment 
consists in retrieving an image sent at the input of the system by measuring the output 
complex field. For focusing experiments A^dof is given by the number of independent input 
segments, which corresponds to the number of degrees of freedom the experimentalist can 
control. Conversely, for detection experiments, the experimentalist cannot access more than 
the number of independent output elements recorded. In the general case, we have A^dof < A^ 
for focusing and A^dof < M. We choose to fix the input (resp. output) pixel size at the input 
(resp. output) coherence length of the system. For this ideal sampling, A^dof is given by the 
input or output dimension of the TM according to the type of experiment : A^dof = A^ for 
focusing and A^dof = M for detection. 

7.2. Experimental Setup 

To measure the TM of an optical system, one needs a way to control the incident light and 
to record the complex optical field. The control is now possible thanks to the emergence 
of spatial light modulators (SLM) and the output complex field can be measured using an 



interferometric method detailed in section 1.3 We want to study an optical system governed 



by multiple scattering in which the propagation function can not be a priori calculated. The 
material chosen is Zinc-Oxyde, a white pigment commonly used for industrial paints. The 
sample is an opaque 80 ± 25 fim thick deposit of a ZnO powder (Sigma- Aldrich 96479) with 
a measured mean free path of 6 ± 2 //m on a standard microscope glass slide. The thickness 
of the sample being one order of magnitude higher than the mean free path, light transport 
inside the sample is in the multiple scattering regime. To generate the incident wavefront, 
the beam from a diode-pumped solid-state single longitudinal mode laser source at 532 nm 
(Laser Quantum Torus) is expanded and spatially modulated by a liquid crystal SLM (Holoeye 
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LC-R 2500). Using a polarizer and an analyzer, we choose an appropriate combination of 
polarizations before and after the SLM to achieve up to a 27r phase modulation depth with 
less than 10% residual amplitude modulation. Only a part of the SLM is modulated and the 



rest of the surface is used to generate a static reference (details in |1.3[ ). The surface of the SLM 
is then imaged on the pupil of a 20x objective with a numerical aperture (NA) of 0.5. Thus a 
pixel of the SLM matches a wave vector at the entrance of the scattering medium. The beam 
is focused at one side of the sample. We image the output pattern 0.3 mm from the opposite 
surface of the sample through a 40x objective (NA = 0.85) onto a 10-bit CCD camera (AVT 
Dolphin F-145B). We do not directly image the surface of the material to mitigate the effects 



of the possible presence of ballistic light. This point is detailed in section [34 

The duration of a typical measurement of the TM is of a few minutes, much less than the 
typical decorrelation time of the medium (more than an hour). From now on, we will ignore 
any time dependence of the medium. It is an important prerequisite for the measurement and 
exploitation of the TM. 



Laser 532nnn 



Reference 
Part 




SLM 



Figure 1. (color online) Schematic of the apparatus. The 532 nm laser is expanded and 
reflected off a SLM. Polarization optics select a phase only modulation mode. The modulated 
beam is focused on the multiple-scattering sample and the output intensity speckle pattern is 
imaged by a CCD-camera. L, lens. P, polarizer. D, diaphragm. 



1.3. Acquisition Method 

Measuring a TM in optics raises several difficulties coming from the impossibility of having 
direct access to the phase of the optical field. Generally, one can have access to the complex 
optical field using interferences, for instance with a known wavefront and a full field phase- 
shifting interferometry method [25 J. For any input vector, if the reference phase is shifted by 
a constant a, the intensity in the m^^ output mode is given by : 
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N 



iZ = l^rT = K + J2^'''k^nK\' (6) 

n 
N / N 



where 5^ is the complex amphtude of the optical field used as reference in the m^^ output 
mode. 

Thus, if we inject the n^^ input mode and measure I^J^^J^ and /m , respectively the 
intensities in the m^^ outgoing mode for a = 0, 7r/2, n and 37r/2 gives : 



/ 7-0 _ T-TT \ I 7^7? — /m 

V-^m -^m; I •_V ^ — ^ h (1\ 

Commercial SLMs cannot control independently phase and amplitude of an incident 
optical field. Our setup uses a phase-only modulator since the phase is the most important 
parameter to control in wavefront shaping f26']. A recent method [i27il allows phase and 
amplitude modulation with similar commercial modulators, nevertheless, we use in this work 
a phase-only configuration to keep the setup as simple as possible. 

Controlling only the phase forbids us to switch off or modulate in amplitude an incident 
mode. Thus we can not directly acquire the TM in the canonical basis. We want to use 
a basis that utilizes every pixel with identical modulus. The Hadamard basis in which all 
elements are either -hi or -1 in amplitude perfectly satisfies our constraints. Another advantage 
is that it maximizes the intensity of the received wavefront and consequently decreases 
the experimental sensibility to noise [28J. For all Hadamard basis vectors, the intensity is 
measured on the basis of the pixels on the CCD camera. We obtain the TM in the input 
canonical basis by a unitary transformation. The observed transmission matrix Kohs is directly 
related to the physical one K of formula [3] by 

Kobs = KX 5ref (8) 

where S^ei is a diagonal matrix of elements where 5^^ = Sm represents the static 
reference wavefront in amplitude and phase. 

Ideally, the reference wavefront should be a plane wave to directly have access to the K 
matrix. In this case, all s^ are identical and Kohs is directly proportional to K. However this 
requires the addition of a reference arm to the setup, and requires an interferometric stability. 
To have the simplest experimental setup and a higher stability, we modulate only 65 % of the 
wavefront going into the scattering sample (this correponds to the square inside the pupil of 
the microscope objective as seen in Figure [T]), the speckle generated by the light coming from 
the 35 % static part being our reference. S^ref results from the transmission of this static part 
and is now unknown and no longer constant along its diagonal. Nevertheless, since S^ei is 
stationary over time, we can measure the response of all input vectors on the m^^ output pixel 
as long as the reference speckle is bright enough. We will show in the next subsection that 



Controlling Light Through Optical Disordered Media : Transmission Matrix Approach 7 

we can go back to the amplitude elements of S^ei and that the only missing information is 
the phase of the reference. We will quantify the effect of the reference speckle and show that 
neither does it impair our ability to focus or image using the TM, nor does it affect our ability 
to retrieve the statistical properties of the TM. 

1.4. Multiple Scattering and Random Matrices 

Under certain conditions, the rectangle A^ by M TM of a system dominated by multiple 
scattering [29] amounts to a random matrix of independent identically distributed entries of 
Gaussian statistics. We will study the properties of such matrices thanks to the Random 
Matrix Theory (RMT) and investigate the hypothesis for which the TM of multiple- scattering 
samples are well described by this formalism. 

We introduced previously that the SVD of a physical matrix is a powerful tool to study 
the energy repartition in the different channels. For matrices of independent elements of 
Gaussian distribution, RMT predicts that the statistical distribution p(A) of the normalized 
singular values follows the Marcenko-Pastur Law OOll . For 7 = M/N, it reads : 



p(A) = -^V(^'-^miJ(ALx-A2) (9) 

zttA 



^min? ^^maxj 



VA G [Amin, An 

With Amin = (1 — a/Vt) the smallest singular value and A^ax = (1 + a/Vt) ^he 
maximum singular value. In particular, for 7 = 1, the Marcenko-Pastur law is known as the 
"quarter circle law" ll3T1l and reads : 



p(A) = -V4-A2 (10) 

Those predictions are valid for a TM of independent elements i.e. without any correlation 
between its elements. This hypothesis can be broken by adding short range correlation but also 
non trivial correlations such as ones appearing when introducing the flux conservation [32J. 
Since we only experimentally record and control a small part of the field on both side of the 
sample, we will not be affected by this conservation condition. 

Another straightforward parameter that could have an influence on the hypothesis of 
independent elements is the size of the segments uses as sources (on the SLM) and receptors 
(on the CCD). If the segment size is less than the coherence length {i.e. the size of a speckle 
grain) two neighboring segments will be strongly correlated. The size of input (resp. output) 
segments is set to fit the input (resp. output) coherence length. This way, the TM K should 
respect the hypothesis of independent elements and follow the Marcenko-Pastur law. 

Nevertheless the method of the acquisition introduces correlation in the TM recorded. 



For M = N, K should respect the Quarter Circle law but the methods detailed in section [L3 
gives only access to the observed TM matrix Kohs = S^^i.K. whose elements read : 
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N 



rnn / j ^mj'^jn '^mn^mm \^^) 

3 

To study the singular value distribution of our TM, we have to be sure that we can go 
back to the precise statistics of K using the observed matrix i^obs- 

The effect of the reference S'ref is static at the output. This means that its influence is the 
same for each input vectors. Calculating the standard deviation of an output elements over the 
input basis vectors we have : 

^«P^ = ^{\kmn?)Js^m\ (12) 

= \l {\kmn?)mn\^rnm\ ^rU G [l, M] 

We made the hypothesis that each input k- vector gives an output speckle of the same 
mean intensity, it reads (|fcmnP)^ = (l^mnP)^^- This is valid in the diffusive regime as 
long as the illumination is homogeneous. We select for the illumination only the center part 
of an expanded gaussian beam. We also neglected the correlation effects that are expected 
to be small regarding our system, such as Memory Effect [l2T1l (see section |4]) and long 
range correlations [331. We can now define a normalized TM Kfu whose elements k^^ are 
normalized by \^{\k^\ 



|2\ 
I m 



yLobs 

kin = .,...."". .„ (13) 



^f¥^. 



mn 



^AlWYr 



Sr. 



(14) 



This matrix is not influenced by the amplitude of the reference and reads 



Where Scf^ is a diagonal matrix of complex elements of norm 1 representing the phase 
contributions of the output reference speckle. In contrast with S^ref, S^j^ is a unitary matrix 
since aS^^M = / where / is the identity matrix. Thus if K = UT^V^ then Kfn = U'T^V^ where 
U' = S(p X U is Si unitary matrix. Therefore, Kfn and K have singular value decompositions 
with the same singular values. 

We show in figure [2] that the statistics of the observed TM Xobs do not follow the quarter 



circle (formula [TO]) law whereas the singular value distribution of the filtered matrix K^\ is 
close to it. Remaining correlation effects are due to inter-element correlations between nearby 
pixels. Taking a sub-matrix of the TM with only one element out of two, its statistics finally 
follows correctly the expected quarter-circle law. Those results are comparable to the ones 
obtained for transmission matrices in acoustics ll29l . It is to be noticed that we do not expect 



absorption to affect the statistical properties of the matrix measured since absorption does 
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not induce any new correlation and since the measure is not sensitive to the total energy 
conservation. 

1^ ' ' 1 
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Figure 2. Singular values distributions of the experimental transmission matrices obtained by 
averaging over 16 realizations of disorder. The solid line is the quarter-circle law predicted 
by RMT. With + the observed matrix, with D the matrix filtered to remove the reference 
amplitude contribution and with o the matrix obtained by filtering and removing neighboring 
elements to eliminate inter-element correlations. 

This clearly demonstrates that we can directly study the properties of K even if we 
only have access to the observed matrix i^obs- We also want to point out that one has 
to be careful when interpreting deviation from a given prediction on SVD results. Once 
the matrix is measured and its singular values analyzed, one has to investigate the causes 
of the deviation from the Marcenko-Pastur law. Those causes could come from physical 
phenomenons involved in light propagation in the medium that break the hypothesis of the 



absence of correlation in the TM (like the apparition of ballistic contributions, see section [34]) 
but could also be inherent to the experimental method. However, Kobs brings all informations 
on singular values that describe the physics of the propagation in a given setup. 

So far, we compared the experimental singular value distributions for a symmetric case 
(M = A^). But what would be the effect of breaking the input / output asymmetric raio M/N 
? In the same setup, we experimentally acquired TMs with a fixed number of controlled pixels 
on the SLM (A^ = 1024) and changed the number of independent output pixels recorded by 
the CCD (M = 7A^) to achieve various ratio values (7 G {1;1.5;2;3;4;5}). Normalized 
singular value distributions of experimental TM filtered are shown in figure |3] and follows 
qualitatively the RMT prediction of equation [TO} 

We see that with increasing 7, the range of A narrows. This effect has been measured in 
experiment in the acoustic field ll34l . Physically this means that when 7 increases, all recorded 
channels tend to converge toward the mean global transmission value A = 1. When M > N, 
we record more independent informations than the rank of the TM and averaging effects lead 
to a distribution that is peaked around a unique transmission value when M :^ N. We notice 
that the experimental distributions deviate from the theoretical ones as 7 increases. To have 
TMs with different values of 7, we record a large TM with 7=11 and then create submatrices 
by taking lines selected randomly among the lines of the original TM. By increasing 7 we 
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Figure 3. Normalized singular value distributions for different asymmetric ratios 7 = 1, 1.5, 
2, 3, 4, 5 and 6. (dash line) experimental results, (solid line) theoretical prediction from RMT. 
Calculations on experimental TM are made using the filtered matrix Kfn. 



increase the probability of having neighboring pixels in the TM and then are more sensitive to 
correlations between nearby pixels. This modifies the statistics and decreases the number of 
independent output segments. This explains the deviation from the theory. Other effects that 



could bring correlations in certain experimental conditions are discussed is section 3.4 



2. Phase Conjugation 

2.7. Phase Conjugation and Degrees of Freedom 

It was demonstrated [[20l |35l that one can take advantage of multiple scattering to focus on 
tight spots, thanks to the reversibility of the wave equation. 

In such an experiment, the responses to short temporal signals emitted by an array of 
sources at a receiver at the output of a disordered medium are recorded. Those signals, 
linked to the Green's functions associated to the couples (source / receiver) are sent reversed 
in time. The waves generated converge naturally towards the targeted spot. The waves are 
focused in space and time. Phase Conjugation is the monochromatic equivalent to the previous 
concept named Time Reversal. Those techniques are robust methods to achieve focusing and 
imaging. In this chapter, we will experimentally and theoretically study the efficiency of phase 
conjugation on focusing. 

When using Phase Conjugation for focusing, it corresponds to shape the input wavefront 
to put in phase the contributions of each input pixel at a desired output target. Denoting E^^^^^ 
the output target vector, the input vector that estimates the desired output pattern £'out^^* using 
Phase Conjugation is given by : 



E""" = KlE'Z^' 



:et 



out 



(16) 



The effective output vector £'°"* will be given by : 



^out ^ j^^in ^ X.K^Ef^^'' 



(17) 
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We will now demonstrate the dependence of the Phase Conjugation accuracy on A^dof 
for the simple case of a single focusing spot. We will compare the efficiency of perfect Phase 
Conjugation (with amplitude modulation) and of phase only Phase Conjugation. The intensity 
output field resulting from the Phase Conjugation focusing on the f^ output pixel reads : 

out|2 






Y.'^mikjil' (18) 

For m ^ j, the average intensity of sum of the the incoherent terms reads : 

IN.N \ 

N \ /N,N 

/ , I ^rnl I Fjd ) + ( Z_^ ^rnl k^i> kji kji> 

= iv(|fc|2)Vo 

For m = j, the average intensity of the sum of the coherent terms reads : 

{\CL/)=((hf^A ) (20) 



N \ /N,N 

'jl' I 



= N{N-l){\k\'f + N{\k\') 
^N^\k\^f ViV>l 

We define the energy signal to noise ratio SNR for a single spot focusing as the intensity 
at the target point over the mean intensity elsewhere after phase conjugation. We directly have 

SNR ^N = A^DOF VA^ > 1 (21) 

Similar results are obtained considering detection instead of focusing. If E'out is the 
output pattern corresponding to an unknown input target mask E^^''^^^ to be detected, similarly 
to equation [171 the reconstruction E'lmage of the target is given by : 



^image = i^^^out = KIK.EI^'^^' (22) 

The reconstructed pattern for detecting a single target on the f^ input pixel reads : 

2 

(23) 



img 



2 



M 



I 
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And the signal to noise ratio for the reconstruction of a single pixel reads : 

SNR ^ M = A^DOF VM > 1 (24) 

2.2. Experimental Phase Conjugation 

To perform such an ideal Phase Conjugation focusing, one has to control the amplitude and 
phase of the incident light. Our experimental setup uses a phase-only SLM to modulate the 
light. Thus the input l^^ component of the input field to display in order to focus on the j^^ 
output pixel reads /c*^/ 1 /c*^ | . The resulting output intensity pattern for a phase only experiment 
reads : 

1..N 7 * l..A^ 

^m I) = y^^rnlTT:^ " 

/ I^J^I I 



1> = E^-^^ = El^-^l (25) 

/ I^J^I I 

Calculating the intensity of incoherent and coherent sums in the same way as 



equations 19 and 20 we have : 



in 



{IsT^jH = N {\k\') (26) 

(ICiiT) -^'(1^1)" V7V>1 (27) 

We define SVRexp the expected signal to noise ratio that reads : 

SVRexp ~ A^DOFTTnlv (28) 

Experimentally, for a setup and a random medium satisfying the hypothesis introduced 



1.4 



(elements of the TM independent with a Gaussian statistic), we have 



IT 



0. 7813611. Finally : 



(IfcP) - 4 



1^^ 



4 



^^Rexp - T^DOF - O.TSTVdof (29) 



one 



To achieve such a 57VRexp using the experimental setup presented in section 1.2 
has to remove, after measuring the TM, the reference pattern. If the reference part is still 
present, the intensity coming from the segments not controlled of the SLM will increase the 
background level. This effect (detailed in Appendix) modifies STVRgxp. In our experimental 
conditions : 

SVi?exp - 0.5iVDOF (30) 



We will now use equations [16] and [22] to experimentally achieve focusing and detection. 
We tested Phase Conjugation for simple targets focusing and for point objects detection 
through the medium. The results are shown in figure |4| 

To experimentally study the effect of A^dof over the quality of focusing, we measure 
the energy signal to noise ratio for a single spot focusing. We represent in figure [5] this 
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Figure 4. Experimental results of focusing and detection using Phase Conjugation. Results 
are obtained for A^ = M = 256. We show typical output intensity pattern for a one spot target 
focusing a. and multiple target focusing b.. We present similar results for the detection of one 
point c. and two points d. before the medium. The point objects are synthesized using the 
SLM. 



experimental ratio as a function of the number of controlled segments A^ on the SLM and 
compare it to the theoretical prediction SNR ^ 0.5A^dof. Those experimental results 
correctly fit the theoretical predictions. 




10 10 10 

Number of independant input pixels 

Figure 5. Focusing Signal to Noise Ratio as function of the number of SLM controlled 
segments N. We show in solid line the theoretical prediction SNR ^ O.SA^dgf and with x 
the experimental results. 



2.3. Comparison with Sequential Optimization 

The first breakthrough using phase-only Spatial Light Modulator for focusing through a 
disordered medium were reported in [8J. The technique was a feedback optimization that 
leads to substantial enhancement of the output intensity at a desired target. We will present 
the principle of this approach and compare it to techniques requiring the measurement of the 
TM. 
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The experimental setup is roughly the same than the one presented in figure [TJ The 
wavefront phase is controlled onto the SLM on A^ groups of pixels, and the full output pattern 
is recorded by the CCD camera. One point target is chosen and the optical field E^^ on this 
pixel is the linear combination of the effects of the A^ segments on the modulator : 

N 
n 

where kmn are the elements of the TM as defined in equation [3] and 0^ the phase imposed 
on the n^^ pixel of the SLM. 

The input amplitude term £"0 is constant for a phase-only spatial light modulator with an 
homogeneous illumination. The idea is, for each pixel, to test a set of phase values and keep 
the one that gives the higher intensity on the target output point. The goal of this sequential 
algorithm is to converge toward a phase conjugation with an efficiency depending on the finite 
number of phase tested [8J. Different phase shifting algorithms were studied and their relative 
efficiencies in presence of time decorrelation are analyzed in [37J, we will not discuss this 
problem in this paper. 

In a single spot focusing experiment for a sufficient number of phase shifts tested, the 
sequential algorithm and the phase conjugation computed by the knowledge of the TM are 
equivalent. Nevertheless, both techniques have different advantages and drawbacks. To be 
able to change the position or the shape of focal pattern, sequential optimization needs a 
new optimization for each new pattern. This forbids real time applications. On the contrary, 
once the TM recorded, all information is known and permits to calculate the input wavefront 
for any output desired pattern, allowing for instance to change the focus location at a video 
rate. Nevertheless, sequential optimization presents two strong advantages. With a setup 
comparable to the one shown on figure [TJ sequential algorithm is more robust to calibration 
errors than the TM measurement since, even if the exact values of the phase shifts tested 
are unknown, the algorithm automatically selects the best phase mask. Another asset of 
this technique is that no interferometry measurement is needed. Since only the intensity 
is recorded, it allows experiments using incoherent phenomenon such as focusing optical 
intensity on a fluorescence probe inside a disordered material ll36ll . 

3. Noise and Reconstruction Operator 

Experimentally, one cannot perform a measurement without noise. In our setup, the main 
noise sources are the laser fluctuation, the CCD readout noise and the residual amplitude 
modulation. Because of the error made on the TM measurement, one has to find, for 
focusing or image detection, reconstructions operators very resilient to experimental noise. 
These operators only give an estimation of the solution. From now on, we will make 
the distinction between the two sources of perturbation on the reconstructed signal : the 
experimental noise inherent to the measurement errors, and the reconstruction noise, inherent 
to the reconstruction operator. We will discuss the influences of those two noises on the 
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performance of the technique used. Due to spatial reciprocity, focusing and image detection 
are equivalent. We will thus concentrate in this section on the case of image detection. 

Inverse problems are widely commons in various fields of physics. Optical image 
transmission through a random system is the exact analog of Multiple-Input Multiple-Output 
(MIMO) information transmission through an unknown complex environment. This has been 
broadly studied in the last decades since the emergence of wireless communications [131 . 
Similar problems have also been studied in optical tomography in incoherent llT6l and 
coherent regime ifTTl . In a more general way, inverse problems has been widely study 
theoretically for more than fifty years. A. Tikhonov ll38l proposed a solution for ill-posed 
system. Tikhonov regularization method was then broadly used in various domains to 
approximate the solutions of inverse problems in the presence of noise. 

3.1. Experimental Noise 

Prospecting an optimal operator for image reconstruction, the first straightforward option is 
the use of inverse matrix K~^^ for a symmetric problem (A^ = M) since i^o^s^obs = I 
where I is the identity matrix. Inversion can be extended when N ^ M with the pseudo- 



inverse matrix 



^L-^obs 



-1 



^obs' i^ol^d i^obg. This operator provides a perfect focusing but 
is unfortunately very unstable in presence of noise. Indeed, the singular values of K~^^ are the 
inverse ones of i^obs- In other words, inverse operators tend to increase the energy sent through 
otherwise low throughput channels. This way, it normalizes the energy coming out from all 
channels. Singular values of Xobs below noise level have the strongest contributions in K~^^ 
but are aberrant. The reconstructed image is hence dominated by noise and bears no similarity 
with the input one. Another standard operator for image reconstruction or communication 
is Time Reversal. This technique has been studied in acoustics ifTSl and electromagnetic 
waves ll39l . Being a matched filter [40J, Time Reversal is known to be stable regarding 
noise level. Unlike inversion, it takes advantage of the strong singular values by forcing 
energy in the channels of maximum transmission. This way, singular values below noise 
level carry few energy and do not degrade reconstruction. Its monochromatic counterpart the 
Phase Conjugation uses the operator Kl-^^ (the operator used in 2.2 for focusing and target 
detection), i^o^gi^obs has a strong diagonal but the rest of it is not null. We will see later that 
this implies that the fidelity of the reconstruction is not perfect and depends on the complexity 
of the image to transmit. 

An intermediate approach is to use a Mean Square Optimized operator. This operator, 
that we will call MSO an denotes W is the Tikhonov regularization Il38l for linear perturbated 
systems mentioned earlier. It takes the experimental noise into account minimizing 
transmission errors, estimated by the expected value E < \W.E'''^^ — E^^'\ \W.E'''^^ — E^^'\ \. 
For an experimental noise of variance No^j on the output pixels, W reads : 



W 



Kl^^-K^i^s + No^.I 



-1 



^o^bs- (32) 



This operator stabilizes inversion through the addition of a constraint depending on the 



Controlling Light Through Optical Disordered Media : Transmission Matrix Approach 16 

noise level. This operator is intermediate between the two previous one. If Noa- = 0, Vl^ 
reduces to the inverse matrix K~^^, which is optimal in this configuration, while for a very 
high noise level it becomes proportional to the transpose conjugate matrix Kl^^, the phase 
conjugation operator. For an intermediate noise level, the channels of transmission much 
greater than the noise level are used 'like in inversion' operator and channels much below 
noise level are used 'like in phase conjugation'. 

3.2. Reconstruction Noise 

The reconstruction noise is the perturbation brought by the reconstruction operator, even 
without measurement noise. This noise is intrinsic to the techniques used. The only operator 
to provide a perfect reconstruction, i.e. without any reconstruction noise, is pseudo-inversion 
(or inversion for M = A^). Nevertheless, this operator is not usable in a general case since it is 
very sensitive to experimental noise. Phase conjugation is stable in presence of experimental 
noise but creates perturbations of the reconstructed signal. As seen in section |2?Tj the signal 
to noise ratio for a single focal is proportional A^dof- The same goes for imaging a single 
spot, the target to detect generates a noise proportional to 1/A^dof on the rest of the image. 
For multiple targets, each target detection adds a perturbation for the reconstruction of the 
others, thus the intensity signal to noise ratio is approximately given by A^f/A^DOP where 
A^f is the number of targets to detect. In other words, the fidelity of the reconstruction with 
Phase Conjugation is not perfect and depends on the complexity of the image to transmit. This 
limitation forbids the reconstruction of a complex image (with A^f ^ N) by Phase Conjugation 
through a symmetric system (with A^ = M) since the reconstruction noise is of the same order 
of magnitude as the useful signal. 

For any reconstruction operator Op, if K is known exactly (not affected by the 
experimental noise), the mean normalized reconstruction noise A^rec for a single target 
detection can be estimated by : 

'''~{\OpxK\U^^^ ^''^ 

Nj-ec is define here by the ratio between the mean intensity {\0p x K\'^^)^^^ brought by 
the reconstruction of a target outside its position and the mean intensity {\0p x K\'^^)^^^ at 
the target position. The average is made over all input target positions. Graphically, A^rec is 
the inverse of the signal to noise ratio defined by the mean of the intensity of the diagonal of 
Op X K over the mean intensity elsewhere. 

To study the reconstruction noise brought by MSO operator for different values of A^o^ we 
show in figure [6] the simulation results for the reconstruction noise as a function of the No^j. It 
is important to understand that, for a given No^j, W is the ideal operator for an experimental 
energy noise level of A^o^. In that study, the simulated TM K is noise free. We see that for 
No (J = {i.e. for inversion operator), the reconstruction noise A^Orec is null and for No^r -^ oc 
{i.e. for Phase Conjugation operator) A^Orec tends to A^dof (here A^dof = M = A^). This 
behavior is obvious when looking at the diagonal of | Vl^ x Kp compared to the others values. 
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Figure 6. Reconstruction noise as a function of No^j. Results are obtained in simulation for a 
TM modelized by a 20 by 20 random matrix of Gaussian statistics. Nocr is normalized by the 
total transmission energy ^^ A^. In inset we show the values of\WxK\'^ with a logarithmic 
scale. 



3.3. Toward Optimal Reconstruction 

In order to detect an input image with the best accuracy, one has to find the good balance 
between the sensitivity to experimental noise and the influence of the reconstruction noise on 
the recovered image. The mathematical optimum operator for those constrains is the MSO 



operator introduced in equation [32| To lower the reconstruction noise, one can only increase 
the parameter A^dof- We investigated two possibilities: averaging, and increasing the number 
of output modes M. 

A possible way to virtually increase A^dof is to perform the transmission of the same 
image through different realizations of disorder. With several realizations, the reconstruction 
noise, depending on the TM, will be different and averaging will lower perturbations on the 
final image. It is the monochromatic equivalent of using broadband signals, which takes 
advantage of temporal degree of freedom Il24l where uncorrected frequencies are used to 
increase A^dof- Iterating with different physical realizations of disorder is restrictive since 
it requires to actively move or change the scattering sample. If we are only interesting in 
detecting an amplitude object, an other way to change the effective TM seen by the object is to 
illuminate it with different phase wavefronts. It is formally equivalent to transmitting the same 
image through different channels as if the image propagated through different realizations of 
disorder. For any amplitude object E^^^ with e^J G [0, 1], the effective incident field on 
the sample is S'^.E^^^ where S'^ is a diagonal matrix containing only phase terms. The output 
complex field pattern for a given TM K reads : 



E^ 



out 



K5;.£;obj = i^^^obj (34) 

With K' = K.S'^. Physically, it is as if the same object E'obj was transmitted through 
a different part of the medium. Another simple way to see it is that by changing S'^ (i.e. 
changing the 'virtual realization') we change the projection of E'obj on the input singular 
modes. Thus the same image propagates through different channels. For a given phase 
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illumination, the amplitude object is estimate by E'img : 



E^ 



img 



\W.E^ 



out I 



\W.K.S'^.Eo^^\ = \W.K\Eo^^\ 



(35) 



With W the optimal MSO operator corresponding to the TM K for a given noise 
level. One now has to change the illumination wavefront S'^ and average the different E'l^g 
obtained to lower the reconstruction noise. To illustrate this effect we show in figure [7] 
simulation results for the reconstruction of a simple amplitude object using two different 
phase illuminations and averaging over 10 different illuminations. The matrix is know exactly 
in the simulation so all perturbations come from reconstruction noise. We see that the two first 
realizations reproduce the amplitude object with different noises with amplitude of the same 
order of magnitude. Averaging over multiple illuminations significantly reduces the noise 
background. 




Figure 7. Comparisons of reconstruction simulations of a simple object for two different 
'virtual realizations' (a. and b.) and the reconstruction averaged over over 10 'virtual 
realizations' (c). Simulations are made for a 100 by 100 random matrix of independent 
elements of Gaussian statistics. The amplitude object consists of three points of amplitude 
one and the rest of the image to zero. 

Since our SLM is a phase only modulator, in order to simulate amplitude objects, we 
generate what we call 'virtual objects'. We use different combinations of random phase masks 
to generate the same virtual amplitude object E^^^ by subtracting two phase objects. From 
any phase mask E^^^^^ one could generate a second mask ^'p^ase where the phase of the m^^ 
pixel is shifted by s'^Ktt. The optical complex amplitude of elements of the second mask are 
linked to those of the first mask by : 



.(2) 



g(l)^giSm^^ 



(36) 



,(i) 



<j) 



with em the nr'^ element of £':{„„„. We have the relation between the two masks and the 



amplitude of the 'virtual object' : 



phase 



|f(2) _fi;(l) I 
I phase phase I 



2. sin (^obj7r/2) 



(37) 



.(2) 



(1) 



it is important to note that the norm of Ei^^^^ — -E^ ^^^ is fixed and controlled whereas 
its phase is random, uncontrolled and differs from a 'virtual realization' to another. The term 
^phase I ^^^ ^^ estimated by : 



-t^amp — l-C'phase 
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Ei^g = |H/.(^S-^S)l (38) 

where E^J^ (resp. E^J^) is the complex amplitude of the output speckle resulting from 

the display of the mask ^'pJase (^^^P- ^phise)- F^^i^ ^ given set of random phase masks, we 



can define in the same way as in equation 34 a diagonal matrix S'i containing only phase 



terms. This term differs from a 'virtual realization' to another. This way, the reconstruct 
image reads : 

^img = |H^.5;.^amp| (39) 

This expression is formally equivalent to equation [35} thus generating 'virtual objects' 
using random phase masks is mathematically equivalent to illuminating a real object by 
random phase-only wavefront. 

3.4. Optimal Reconstruction in Presence of Ballistic Light 

We experimentally used this technique to virtually increase A^dof, and we average the results 
to lower the reconstruction noise. The complex amplitude image to detect is a grayscale 32 
by 32 pixels flower image. The system is symmetric with N = M = 1024. The first step to 
reconstruct the input image with the best accuracy is to find the variance of the experimental 
noise No^r in order to obtain the optimal MSO operator. To that end, we numerically 
iterate the reconstruction operation for different values of the constrain in the MSO operator 
(equation [32]) and keep the one maximizing correlation between the reconstructed image and 
the image sent, hence obtaining an estimation of the experimental noise level. This step is 
made a priori knowing the image sent but could be performed by any technique leading to 
an accurate estimation of the noise level. We detailed in filSil the results in a purely diffuse 
system where ballistic contributions have no measurable effect. We study here the results in 
presence of quantitative ballistic contributions. 

In previous experiments, the focal planes of the two objectives are not the same. If 
the two focal planes coincide, if we remove the sample, an incident k- vector matches an 
output k- vector. Thus a pixel of the SLM corresponds to a pixel of the CCD. In the same 
configuration with a scattering medium, if the part of the energy ballistically transmitted is 
important enough, this baUistic part can be isolated in the TM. 

To observe those contributions and study the accuracy of the reconstruction methods 
in that particular case, we experimentally use a thinner and less homogeneous area of 
the medium with a symmetric setup (M = A^) where input and output segments are 



approximately the same size. We test the different reconstruction methods presented in [33 
with one realization and averaging over 40 'virtual realizations' with random phase masks. 
Results are shown in figure [8j As expected, we see that the inversion operation does not allow 
image reconstruction, even with averaging, whereas optimal MSO allows a 93% correlation 
for 40 averaging (and a 60% correlation for one realization). Phase conjugation which is 
supposed to be stable in presence of noise gives a weak correlation of 31% with averaging. 
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Figure 8. Comparisons of the reconstruction methods in presence of balHstic contributions. 
a. initial gray scale object and b. a typical output speckle figure after the opaque medium, c, 
e. and g. are experimental images obtained with one realization using respectively inversion, 
phase conjugation and MSO operator, d., f. and h. are experimental images averaging over 
40 'virtual realizations' using respectively inversion, phase conjugation and MSO operator. 
Values in insets are the correlation with the object a.. 



Those results are similar to the ones obtained when we measure the TM of a system 
governed by multiple scattering ifTSll for inversion (11% correlation with averaging) and for 
MSO operator (94% correlation with averaging) but are very different for phase conjugation 
(76% correlation with averaging). We see in figure [8]^e.) and (f.) that most of the energy of 
the reconstructed field is localized on very few points. In presence of strong enough ballistic 
components, the direct paths associated give rise to channels of high transmission linking one 
input k- vector to the same output k- vector. Figure [9]^a.) shows normalized singular values for 
both experiments and for a simulated random matrix following the quarter circle law. For the 
'ballistic' experiment, two high singular values can not be explained with a TM governed by 
multiple scattering only. We show for comparison in figure [9]^b.) the spatial aspect of the input 
singular vectors corresponding to the first and the 100^^ highest singular values, which should 
not be affected by balistic light. The 100^^ singular vector is spatially randomly distributed in 
energy (as expected for random multiple scattering) whereas the first singular vector mostly 
contains energy in one k-vector. This confirms that, in this particular case, the channel of 
maximum transmission is mostly dominated by the ballistic contribution. 
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Figure 9. (a.) normalized singular values A sorted by decreasing order for a thin non 
homogeneous region of the sample with experimental conditions sensitive to ballistic waves 
(dash line) and for a thicker region of the medium with the two focal planes objectives 
separated by 0.3 mm (solid line). With • we show the normalized singular values of a simulated 
matrix of independent Gaussian elements following the 'quarter circle law'.(b.) amplitude 
of the singular vector corresponding to the maximum singular value (left) and to the 100^^ 
singular value (right). Both share the same amplitude scale. 



Phase conjugation maximizes the energy transmission in channels of maximum 
transmission [40J. Therefore, baUistic high singular values contributions will be predominant 
in phase conjugation, regardless of the image E'obj. Since the ballistic singular vectors are 
not homogeneously spatially distributed in energy, they will not efficiently contribute to an 
arbitrary image reconstruction. MSO is not affected since lower the weight of channels which 
do not efficiently contribute to the image reconstruction, including the ballistic ones. 



3.5. Pseudo-Inversion 

We showed in section |1.1| that when we set the output pixel size to fit the output coherence 
length, the number of degrees of freedom for detection is equal to the number M of segments 
read on the CCD. This way a second approach to increase A^dof and thus the quality of the 
reconstructed image is to enlarge the output observation window. An important advantage is 
that the limiting time in our experiment is proportional to the number of input basis vectors 
A^, thus, in contrast with a focusing experiment where A^dof is equal [8J to A^, increasing 
the degrees of freedom for image reconstruction does not increase the measurement time. 
Another consequence of increasing the size of the image recorded, that is, of increasing the 
asymmetric ratio 7 = M/N > 1, is that the statistical properties of the matrix are deeply 
modified. We saw in section [O] that increasing 7 the range of the normalized singular values 
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decreases OOllHl. A direct consequence is that the smallest non-zero singular value increases 
and reads A^ = (1 — ^1/7). Controlling less information on the input than measured on the 
output results in an averaging effect. An interesting consequence is that for a reasonable 
experimental noise level, one can find a value of 7 for which all singular values (and thus all 
channels energy transmission) are above this noise level. In such a situation, the TM recorded 
is barely sensitive to the experimental noise and since no singular values are drowned in the 
noise, pseudo-inverse operator can be efficiently used. 

To highlight this effect, we experimentally measure the TM of our sample for different 
values of 7 > 1 . For each TM we tested a single realization of the reconstruction of a complex 
image using optimal MSO operator and pseudo-inverse operator. The results are shown in 
figure [TOl^a.). For 7 = 1 we are in agreement with previous results since optimal MSO gives a 
partially accurate image when results of inversion are totally drown in noise. Increasing 7 and 
thus A^DOF strongly improves the quality of the reconstructed image for both techniques and 
reaches a > 85% fidelity for the largest value of 7 = 11, without any averaging. As expected, 
the curves of fidelity reconstruction as function of 7 for both techniques become identical 



when the value of A^ reaches noise level (figure 10 ^b.)). Above this value, pseudo-inversion 
and optimal MSO are totally equivalent. One notice that experimental A^ are always smaller 
than their theoretical predictions. We explain this deviation by the correlation induced in the 



matrix by the reference pattern I^Sref | (see section 1.4). 

We saw in this part that increasing the degrees of freedom and using an appropriate 
operator that takes noise into account permits to drastically enhance the reconstruction 
accuracy of an image detection. Similar results can be obtain for focusing or for output 
wavefront shaping with the same techniques. Such experiments are more difficult to carry 
out since they require amplitude and phase modulation. Increasing 7 = N/M and thus A^ in 
that case also presents the drawback to significantly increase the measurement time. 

4. Memory Effect 

To define the statistical properties of the TM of multiple scattering media in part |1.4| we did 
not take into account weak correlation effects that could break the hypothesis of independent 
elements of the TM. The strongest (first order) of those effects that can be recorded for system 
dominated by multiple- scattering (no ballistic contribution) in the diffusive regime is the 
so-called Memory Effect (ME). It brings correlation with characteristic length L^e which 
depends, in transmission, on the geometry of the observation system EUl . If the speckle grain 
size I/gpeckie is of the same order of magnitude or smaller than the ME range L^e, ME brings 
additional correlations in the TM. 

We define Lspeckie as the full width at half maximum of the autocorrelation function of a 
speckle intensity pattern and L^e as the attenuation length of the ME. Both length are defined 
for a given observation plane. 

The speckle size Lgpeckie depends on the distance between the back of the sample and the 
observation plane D, the width of the illumination area W and the wavenumber k. It reads : 
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Figure 10. Influence of the number of output detection modes, a. correlation coefficient 
between E^^g and £^obj as a function of the asymmetric ratio 7 = M/N for MSO (dashed 
fine) and for pseudo-inversion (soHd line), without any averaging. Error bars correspond to 
the dispersion of the results over 10 realizations, b. Experimental (solid line) and Marcenko- 
Pastur 1301 predictions (dashed line) for the minimum normalized singular value as a function 
of 7. The horizontal line shows the experimental noise level No^^^. 
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In contrast, L^e depends on k, D and the thickness e of the sample and reads : 

ke 



(40) 



(41) 



Therefore, for Vl^ :» e we have L^e ^ ^speckle- This kind of geometry can be convenient 
since it allows to move a focal spot along a distance greater that Lgpeckie by just tilting the 
SLM[I4T1. 

We experimentally fixed D = 1.5 ± 0.2 mm, W = 0.5 ± 0.1 mm and the sample has a 
thickness e = 80 ± 25 //m. For those conditions, we expect Lgpeckie ^ ^me ^1.5 /xm. The 
TM is measured. Each column of the TM is the complex amplitude response of a plane wave 
with a given incident angle. One can easily reconstruct the intensity patterns for different 
incident angles and calculate the correlation function between two patterns : 
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1^71 



k 



{M-m)n\ 



n—n'),j 



zZm \^mn\ X Z/m F^ 



(42) 



{n—n') 



Where ®^ denotes the spatial convolution over the variable m and j is the variable of the 
correlation function corresponding to a displacement in pixel in the focal plane. We average 
results over all combinations of n, n' with n — n' = constant. Results are shown in figure 1 1 
It is clear that the output speckle is coherently translated for distances greater than I/gpeckie- 
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Figure 11. Spatial correlation functions of output speckle for a variation of the input incident 
angle. The correlations functions are drawn as function of the spatial coordinate corresponding 
to the displacement of the speckle pattern in the observation plane. Arrows show the maximum 
of correlation functions for the highest displacements. Those peaks are barely observable due 
to an exponential decay of the maximum correlation value. 



^speckle is deduced from the mean autocorrelation function of the speckle pattern 
f{n-n^)=o,j as response to an incident wave plane. We found Lgpeckie = 1.3 ± 0.48 /xm. We 
define C{x) the correlation coefficient (maximum of the correlation function) between two 
output speckles shifted by x in the observation plane. Exponential decay of the correlation 



C{x) 



oc e 



-x/ Ln 



is generally observed for an open geometry ll2T1l . The measured C{x) 



(figure [TT] and [T2|) fit an exponential decay with L^e = 1-72 ± 0.13 fim. 

We found the expected correlation lengths with a good order of magnitude for autocorre- 
lation (speckle grain size) and memory effect and we were able to qualitatively separate both 
effects. 



Conclusion 



We reinterpreted in this article previous works on optical phase conjugation in term of linear 
system defined by a monochromatic transmission matrix. Within this formalism, we show 
how to take advantage of multiple scattering to achieve image transmission and wavefront 
shaping. We also pointed out how to improve the quality of measurements in the presence of 
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Figure 12. Spatial correlation decay of the memory effect, (solid line) experimental curve 
of correlation coefficient as function of the spatial displacement, (dash line) exponential fit 

C{x) = e-^/^-- with Lme = 1-72 fim. 



noise, taking benefit of Random Matrix Theory and information theory. We finally showed 
how to take benefit of the statistical properties of the TM to go back to physical phenomenons. 

Acknowledgments 

This work was made possible by financial support from "Direction Generale de 1' Armement" 
(DGA), Agence Nationale de la Recherche via grant ANR 2010 JCJC ROCOCO, Programme 
Emergence 2010 from the City of Paris and BQR from Universite Pierre et Marie Curie and 
ESPCI. 



Appendix 

To take into account the effect of the reference beam in the effective SNR, we define F < 1 the 
ratio of the intensity coming from the controlled part over the total intensity. This way, F is 
also the ratio of the area of the controlled part of the SLM over the total illuminate area. One 
way to understand this situation is that A^tot independent segments of the SLM are illuminated 
and we control only A^ = F. Atot of them. We injected the ratio F in equation [19] and [20j SNR 
now reads : 



SNR 



7r/4A2 + (Atot - A) 



TV 



/m^ 



Ntot 



r) 



N/T 



4 



(A.1) 



VA^>1 



For the setup shown in figure [TJ the illuminated area on the SLM represents the 
circumcircle of the square part that we control. The segments illuminated but not controlled 
are used to produce the reference pattern. For this geometry, with an illumination disk of 
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radius R we have F = ^ ^o — = -• The effective SNR corrected to take into account the 
reference beam is therefore : 

SNR ^ Q) ^TVdof = 0.5iVDOF (A.2) 
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